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Abstract 

The effect of nonlinear feedback control strategies on the planform of convection 
in a Boussinesq fluid heated from below is investigated. In the absence of the 
control, given that non-Boussinesq effects may be neglected, it is well known that 
convection begins in the form of a supercritical bifurcation to rolls. Non-Boussinesq 
behaviour destroys the symmetry of the basic state, and through a subcritical 
bifurcation leads to the formation of hexagonal cells. Here we discuss the influence 
of regulation of the lower surface temperature by means of a control mechanism, 
made up of a combination of a proportional linear and nonlinear controller, on the 
stability of the hexagonal cell pattern. 
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1 Introduction 


Our concern is with the effects of feedback control on the evolution of convective insta- 
bilities in a Boussinesq fluid heated from below. The stability of Benard convection is of 
both mathematical and physical importance, and is of particular relevance to material 
processing applications, where it is often more desirable to operate at Rayleigh numbers 
in excess of the critical value above which thermal convection normally occurs, but in the 
absence of such convective instabilities. 

Currently methods for suppressing convection include, damping fluid motion in elec- 
trically conducting materials with magnetic fields, and also operating in a reduced gravity 
environment in order to minimise the effects of buoyant convective currents. 

It is well known, that for non-Boussinesq fluids, where variation of temperature with 
viscosity is taken into account, the onset of convection in a steady fluid layer, heated 
from below, takes the form of hexagonal cells with the flow direction at the cell centre 
determined by the sign of the change in sign in viscosity variation. (See Segel & Stuart 
(1962), Palm et al. (1967), Segel (1965), and Busse (1967).) This cell pattern is induced 
by a subcritical bifurcation to convective rolls. Further increase in the Rayleigh number 
cause a destabilisation of the hexagonal pattern, and roll modes are established. However, 
if the Rayleigh number is subsequently decreased to a value below that at which the roll 
modes first appeared, hysteresis effects are seen, that is, the hexagonal cells are not re- 
established until a value below that at which transition to the roll mode first occurred. 
This is a consequence of the subcritical bifurcation, and an additional consequence due 
to the nature of this bifurcation is that the onset of convection occurs discontinuously 
when the Rayleigh number is increased. 

Considerable research has been done on stratagies to delay the onset of Benard con- 
vection in a fluid layer heated from below and/or cooled from above through modula- 
tion of the boundary conditions and, in particular, the time-periodic modulation of the 
temperature difference across the layer - see Davis (1976), and Donnelly (1990). Pre- 
vious authors have found that this leads to a small increase in the critical Rayleigh 
number (see Meyer, Cannell and Ahlers (1992)). However, Roppo, Davis & Rosenblat 
(1984), found that periodic modulation of the lower surface temperature induces symme- 
try breaking effects sufficient to cause a subcritical bifurcation and thus the no-motion 
state may only be stable for small perturbations. 

Recently, Kelly and Hu (1993), showed that a time-periodic non-planar shear flow sig- 
nificantly stabilises the Benard problem. In a subsequent paper Hall and Kelly (1995), 
it was found that steady and unsteady shear flows removed the preference for hexagonal 
cells at the onset of convection in a non-Boussinesq fluid. This was achieved by suffi- 
ciently splitting apart the critical Rayleigh numbers for each member of a triad of roll 
disturbances which interact to produce the hexagonal cells. 

The mechanism found by Kelly and Hu (1993), is possibly explained by results for 
convection in unidirectional steady flows, see the review by Kelly (1993), and this possi- 
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bility is summarised in Hall (1995). Kelly and Hu found that the maximum stabilising 
effects were for low frequencies. For this reason, Hall (1995) concentrated on the low 
frequency regime, and investigated the interactions between all possible roll disturbances 
for a given value of the Rayleigh number. Note that for small amplitude unsteady flows, 
subcritical convection of any form is not possible, Hall and Kelly (1995), which is in con- 
trast to the results of Roppo, Davis and Rosenblat (1984). Hall (1995), found that there 
exists a threshold frequency, below which imperfections in the system play a crucial role, 
causing the system to respond in a quasi-steady manor. The non-linear problem becomes 
singular, in the sense that, the solution to the problem differs by an 0(1) amount from 
the imperfection free solution. However, even for frequencies small enough to induce 
quasi-steady behaviour within the fluid layer, at any instant in time the convection cells 
do not line up with the direction identified as most unstable on the basis of linear theory, 
that is, the behaviour for such frequencies cannot be predicted from the corresponding 
steady case. The convection patterns in this regime are found to change with time, and 
are extremely complicated with straight rolls existing only for part of a period. It was also 
found that the Kelly-Hu mechanism persists into the nonlinear regime until imperfections 
start to play a role in the system. Thus practical application of such a stabilising method 
must ensure that the frequency of the oscillations required for a full stabilising effect is 
shorter than that over which imperfections have an effect on convection. In addition, the 
necessary horizontal motions may also have an adverse effect in materials processing. 

In various experimental and theoretical investigations, by Singer, Wang and Bau, (see 
references), the bifurcation structure associated with convection has been successfully 
changed by means of a feedback controller in a thermal convection loop, heated from 
below and cooled from above. Thus the Rayleigh number above which convective motion 
occurs was higher than that for the uncontrolled system. In addition, they also induced 
chaotic behaviour for flows which would otherwise be laminar. 

Tang and Bau (1993), implemented similar ideas when investigating the Lapwood 
problem for a horizontal saturated porous layer, confined in a box heated from below, 
and cooled from above. Lapwood convection is relevant to transport processes in the 
mushy region of solidification, and also in thermal porous insulators. The convection in 
such a medium, may be modelled for low to moderate Darcy-Rayleigh numbers by the 
Darcy-Oberbeck-Boussinesq equations. They concluded that the no-motion state could 
be maintained by means of feedback control strategies, for Rayleigh numbers exceeding 
the critical values for the uncontrolled system. In addition, the control mechanism was 
found to leave unaltered the no-motion solution of the classical equations, but changes 
the stability characteristics by altering the system dynamics. The control mechanism 
implemented was either of the proportional, or differential kind or a combination of the 
two. 

The problem under investigation here, concerns the effects of a nonlinear feedback 
control mechanism, that is we modulate the lower surface thermal boundary conditions 
with reference to the midlayer temperature with a combination of proportional and non- 
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linear controllers, the latter being dependent on the square of the reference temperature. 

The aim of this study is to determine firstly if subcritical convection can take place, 
and investigate the stability when nonlinear effects are taken into account. 

In the first instance, we shall derive a system of amplitude equations, which are 
generalisations of those found by Segel (1965), for the problem under consideration. In 
section 3, we will investigate the bifurcation structure of these equations. And finally, in 
section 4 we shall draw some conclusions. 

2 Formulation 

Consider a fluid occupying the region 0 < z < d, and of infinite horizontal extent. The 
fluid is Newtonian, and we apply the Boussinesq approximation. The wall z = 0 is 
maintained at a temperature T 0 + AT, with the upper wall at a temperature T 0 . We 
introduce non-dimensional variables, scaling length, time, velocity, pressure and temper- 
ature according to d, d 2 /c, (agKATdv~ 1 )* , (pQ/cagi/ATd -1 )^ , and AT respectively. The 
governing equations are then as follows, 


v f + i?2v-Vvj = V 2 v — Vp + R 2 Tk, 

(2.1a) 

o' 

II 

> 

> 

(2.1b) 

T t + r?v • vr = v 2 r, 

(2.1c) 


where, a — vf k, is the Prandtl number and R = agATd 3 j kv is the Rayleigh number. In 
the above, subscripts denote partial differentiation, and the operator V 2 is defined as, 

2 d 2 d 2 d 2 

dx 2 dy 2 dz 2 

The velocity, pressure and temperature fields are denoted by the quantities, v = (u, u, u;) , 
P and T respectively, whilst k, represents the unit vector in the positive z direction. Also, 
p 0l u^K and a, are typical values for the density, kinematic viscosity, thermal diffusivity 
and coefficient of thermal expansion respectively. 

Before proceeding, we first describe the feedback control mechanism introduced by 
Tang and Bau (1993). It consists of sensors placed at some horizontal cross-section, in 
our case, the midlayer, and actuators on the lower surface, which react to any deviation of 
the fluid temperature from the desired conductive values by modifying the temperature 
distribution there, and thus, increasing the dissipation of these disturbances. 

The no-motion state of the classical problem is a solution of the equations of motion 
for all values of the Rayleigh number . However, as a result of imperfections within the 
system, the temperature field may deviate from its conductive value, and thus buoyant 
forces are produced, inducing fluid motion. Below the critical Rayleigh number for the 
classical problem, thermal conduction is sufficient to dissipate these disturbances and 
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hence, the no-motion state is restored. The controller introduced, may assist in this 
dissipation, and thus the critical value of the Rayleigh number maybe be increased. 

We assume that the temperature distribution at 2 = is known, or maybe calculated 
to sufficient accuracy and is a continuous function of (x,j/). In the no-motion state, 

T z= o = 0. 

Given some deviation from the initial distribution, the actuators modify the lower surface 
temperature according to, 

T z=0 = k p T+K n T 2 , (2.2) 

where T = T(z = |) and k p and k n are the gains of proportional and nonlinear con- 
trollers respectively. 

Once the disturbance has been dissipated, T z= o again takes its initial value of zero. In 
normal operating conditions, the controller reacts before the disturbance has a chance to 
grow, and any deviation in the lower surface temperature will be small. Thus we make 
the assumption that the control parameters, K p and K n are small, and thus introduce, 

K p = K p 8 , K n = K n 8 , 

where, 6 is a small parameter, and K p , and K n are both 0(1) quantities. 

The walls at z = 0 and z = l are taken to be planar, stress-free surfaces, and con- 
ditions of no slip and zero stress are imposed, together with the assumption that the 
boundaries are isothermal. These conditions may be re-written in terms of the tempera- 
ture variable as. 

On z = 0, 

T = K p T8+K n T 2 8 , 

V 2 T = V 4 T = 0 

and on z = 1, 

T = V 2 T = V 4 T = 0 

Again, T, is defined as the temperature on the surface, z = 

We note that boundary conditions for T are, 

D 2 T = D 4 T = 0, on 2 = 0, 1 

together with suitable conditions on T at the upper and lower surfaces. However, due to 
the control mechanism under investigation, these conditions depend on 8. In the above, 
D denotes partial differentiation with respect to z. 

We wish to concentrate on the nonlinear evolution of interacting modes, and we seek 
solutions with R differing from its critical value by 0(e 2 ), where e is considered to be 
small. We therefore expand R as, 

R = Ro + c 2 f?2 + . . . . (2.5) 


(2.3a) 

(2.3b) 

(2.4) 
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Note that in the above, Ro and i2 2 , are dependent on 6, and in addition, we scale the 
time dependence such that t = e 2 r, where r is an 0(1) quantity. 

Here, we wish to concentrate on the roll of feedback control on the pattern selec- 
tion problem in the weakly nonlinear regime, and apply the results of previous authors 
(see Segel and Stuart 1962, Segel 1965) for the slightly non-Boussinesq fluid to aid in 
interpreting our results. The latter authors considered pattern selection involving rolls 
and hexagons, and found that the problem reduces to the investigation of the non-linear 
interactions of a triad of modes with wave-numbers, 


( k k ) ~ -(- — 

[ x, y) ~ y/2 \ 2 ’ 2 

(k k) - -(- — 
1 x ’ y) ~ y/2\r 2 

{k x ,k y ) = ^( 1 , 0 ). 


(2.6a) 

(2.6b) 

(2.6c) 


In the absence of the controller, we expand the temperature perturbation as, 


T = 1 — z + e sin xz ( XE\ + YE 2 + ZEz) + c.c. + . . . , (2-7) 


where, 



(2.8a) 


(2.8b) 

(2.8c) 


and X , Y and Z, are functions of time only whilst c.c. denotes the complex conjugate. 

If the non-Boussinesq effects are taken to be small, O(e), then the amplitude equations 
for X, Y and Z, are of a similar form to those of Segel and Stuart ^96^, where variation 
of the fluid properties with temperature is taken into account. However, since we assume 
that viscosity is constant, the quadratic term in these equations vanish. We leave the 
discussion of this system of equations until the following section. 

In order to determine the effects of feedback control upon pattern selection develop- 
ment, we shall expand the temperature and velocity fields as follows, 


T = Te + eT^ + e 2 T 2 + ..., (2.9a) 

v = eV] + e 2 v 2 + (2.9b) 


We now substitute the above into the governing equations (2.1a-c) and equating like 
powers of e and then S we obtain a solution in the following manner. 
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In the basic state, the fluid is at rest. The temperature field, Tb ( z ), and Rayleigh 
number Rq are expanded below, 



Tb — Tbo + $Tb\ 4- & 2 Tb2 + - - * , 

(2.10a) 


Ro = Roo + $Rcn + . . . , 

(2.10b) 

where, 

£ 

o 

ii 

H— * 

1 

(2.11) 

and 

Tbi = 7i(l -z). 

(2.12a) 


7i = KpT bo + KnT B0 . 

(2.12b) 

Similarly, 

T B 2 = 72 ( 1 - 2 ), 

(2.13a) 


72 = K P T bi + 2K n T bqT bi- 

(2.13b) 

At 0(e), 

we again expand the temperature and velocity perturbations 

in powers of 6, 


Ti = Tjo + STu + . . . 

(2.14a) 


Vi = V 10 + £v n + .. . 

(2.14b) 

The problem for Tjo reduces to, 



- V 6 T lo + /?ooV?T 1 o = 0, 

(2.15) 

where 

V72 & d 2 

V; = 4 - 

1 dx 2 dy 2 ' 


with the boundary conditions: 



T\o (o) = d 2 t w {0) = d 4 t 10 (o) = o, 

(2.16a) 


iio(i) = r> 2 r 10 (i) = D 4 r J0 (i) = 0 . 

(2.16b) 

The solution to the above maybe written as, 



T 10 = <f)S\n 1 TZ 

(2.17) 

where, 

vfr+Y* = °» 

(2.18a) 


<t> = X{t)Ex + Y{t)E 2 + Z{t)Ez + c.c., 

(2.18b) 
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where, E U E 2 and E 3 are as previously defined in (2.8a-c). In this way we find the critical 
Rayleigh number , Roo to be equal to 277 t 4 /4. 

Again at 0(e 2 ), we expand T 2 and V 2 , as power series in 6, as in (2.14a,b), and find 
that the leading order problem for T 2 q and v 2 o has a solution given by, 



X*20 — ^ ^ Lm^m sin 2>7TZ) 

ra=0,4 

W 20 = Pm^mSm27TZ, 

771 — 0 A 

(2.19a) 


(2.19b) 


(“ 20 ) = £ «-(!“) «»2«. 

Vv 20 / m= 0,4 \®rnyj 

(2.19c) 

In the above, L m ,P m and are constants, determined from the system of equations 

governing this order and the subscripts, x and y denotes partial differentiation with 
respect to the x and y co-ordinates respectively. The functions $ m are given by, 

$0 

= 2(\X\ 2 + \Y \ 2 + \Z\ 2 ), 

(2.20a) 


= 2 YZEi + 2 ~XZE 2 + 2XYE 3 + c.c, 

(2.20b) 

$ 2 

= 0 , 

(2.20c) 

$3 

= 2XYE l E 2 + 2XZE 1 E 3 + 2YZE 2 E 3 + c.c, 

(2.20d) 

$4 

= X 2 E 2 + Y 2 E 2 2 +Z 2 E 2 + c.c , 

(2.20e) 


and thus, 

7T 2 

V 2 4> m = 

We note that this solution satisfies homogeneous boundary conditions on the temperature 
at the lower surface. 

At 0(eS ), we find that Tn satisfies, 


-V*T„ 

+ RotVjT'n = — (R0071 + R01) VjTio, 

(2.21a) 

r„(i) = v ! r„( 0 ) 

r„(0). = (K p + K n )<f>, 

= v 2 r„(i) = v 4 r n (o) = v 4 r„(i) = o, 

(2.21b) 

(2.20d) 

with solution, 

Tn = 

< if{z ), 

(2.21a) 

T(z) = 

At (z — 1) costtz + ^ exp XjZ, 

j= 1,4 

(2.21b) 


where At is a constant given by, 



and Xj are roots to the equation, 



other than ±Z7T. Values for At, and Bj in (2.21a) are found from the boundary conditions, 
which reduce to solving a matrix equation of the form, 

Px = (K p + K n ) Q. (2.23) 

Hence, from (2.22), a value for Rq\ may also be determined. 

Since we wish to consider the evolution of the amplitudes, X, Y and Z, at 0(e 3 ), and 
include nonlinear effects from both the controller and terms in the governing equations, 
we require that e and S are of comparable orders of magnitude. Thus, we shall combine 
the solutions to the problems at 0(e 3 ), 0(i 2 8) and 0(t8 2 ). In other words we now set 

e = 6. 

At this order, the solvability condition required for a solution of the differential system 
to exist yields the following equations for the amplitudes X, Y and Z, 

EX = (e + a) X — aYZ — X (f?i | X | 2 +P | Y | 2 +P \ Z | 2 ) (2.24a) 

EV = (t + a)Y-aXZ-Y(R,\Y\ 2 +P\X\ 2 +P\Z\ 2 ) (2.24b) 

EZ = (e + a)Z-aXY~Z(R 1 | Z | 2 +P | X | 2 +P | V | 2 ) (2.24c) 

In the above, a and a are functions of the control parameters K v and K n and given 

by, 


a - + Y+^j 

+ y (7i A>. + 7 2 «oo) + ( K„t» (i) + K„ (f u (i) + 7l 

+ 24f-3T 2 r(0)-^(-^+ ^ 

L 16 Ri 0 \ 37r ,=1,4 ) 

- 97 r 5 X n . 


(2.25a) 


(2.25b) 


We note that from, (2.23), the constants, Aj, and Bj, are proportional to K p -\- K n . From 
this, and the above expression for a, we obtain the following, simplified expression for 
the dependence of a on the control parameters, K p and K n , 

a = a-i (K p + K n ) + a 2 K n , (2.26) 
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where a x and a 2 are constants. 
In addition, 

and, 

and /?] and P are constants. 


3 Solution of the amplitude equation 

We seek steady solutions of (2.24a-c) with, 

X — x 0 (3.1a) 

Y = y o (3.1b) 

Z = z 0 . (3.1c) 

In the above, x 0 ,yo and z 0 , are constants satisfying the following system of equations, 

0 = (e + a) x 0 - ay^z 0 - x 0 (Ri | x 0 | 2 +P | y 0 | 2 -\-P \ z 0 | 2 ) (3.2a) 

0 = (e + a)y 0 - ax^z 0 - y 0 | y 0 | 2 +P \ x 0 | 2 +P \ z 0 | 2 ) (3.2b) 

0 = (e + a) z 0 - ax 0 y 0 - z 0 | z 0 | 2 +P \ x 0 \ 2 +P \ y 0 | 2 ) (3.2c) 


With the control switched off, K p = K n = 0, we note that a and a are both zero, and 
from Segel and Stuart (1962), and Segel (1965), possible equilibrium solutions are, 


I Conduction : Xo = yo = zo = 0, (3.3a) 

IIa,b Rolls : x 0 = yo = 0, z 0 = ± (^-) 2 , (3.3b) 

IIIa,b Hexagons : x 0 = t/o = z 0 ,z Q = qF (2 T) _1 \/4eT, (3.3c) 

IVa,b Hexagons : x 0 = yo = —zo, zo = T ( 2T )~ 1 \/AeT , (3.3d) 

V Mixed : z 0 — 0, x 0 = yo = ±VRe, (3.3e) 

where Q — P — R\, 4 R = P + R\ and T = P + 4 R. 


The stability of the different solutions may be found by an examination of the lin- 
ear growth rates of small perturbations to x 0 , y 0 and zq as discussed in appendix 4 of 
Segel (1965). The equilibrium forms for zq given by (3.3a,e), are shown in figure 1, with 
the Prandtl number set to unity. Unstable solutions are represented by dashed curves. 



— — ( R20 + Ru + R02) ■ 
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Figure 1 : Equilibrium forms for z 0 , given by (3.3a,e) 


We now go on to describe how the bifurcation picture is modified by the effects of 
feedback control on the system. We first write down the equilibrium solutions satisfied 
by the controlled problem, 


I Conduction 

H 

o 

II 

o 

II 

£ 

II 

O 



(3.4a) 

IIa,b Rolls 

•• zo = t/o = 0, z 0 = ± ) 

/ 

1 

2 

i , 


(3.4b) 

IIIa,b Hexagons 

: xo = yo = zo, zq = (2 T)~ x (- 

-a TV 

/ 

^a 2 + AT (e + a) \ , 

(3.4c) 

IVa,b Hexagons 

■ xo = yo = - Z 0 , Zo = (2T)~ 1 1 

[-a=F 

\Ja 2 + AT (e + a) 

) i (3.4d) 

V Mixed 

—a i ( 

: z o — ~~q ? = 2/o — i /? 2 1 f + o — 

Ria 2 \> 

Q 2 ) ’ 

(3.4e) 
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where Q , P,R\,R and T are as previously defined. 

Clearly, for solutions of types II-V to exist, t must be greater than some threshold 
value given by, 


III, IV 

a 

: e = Bl = — — a, 
4 T 

(3.5a) 

II 

: t = B2 = — or, 

(3.5b) 


R i<z 2 


V 

: .-Al= ^ 

(3.5c) 


In addition, it can be shown by a linear stability analysis, that the hexagonal cell pattern 
loses stability at, 

t=A2= ^£- a - (x6) 

The first point to note is that a subcritical bifurcation to the hexagonal cell pattern 
remains viable. This is due to the fact that all three modes will have the same critical 
values of e, corresponding to linear instability, since the control mechanism has no pre- 
ferred direction, unlike say the shear flow mechanisms discussed by Hall and Kelly, where 
the effect of the shear is felt only by modes not aligned parallel to the shear. 

We now go on to examine the case where only the proportional controller, K pi is in 
operation, and see how the bifurcation diagram changes. 

Figure 2 illustrates how the bifurcation picture is modified for K p = ±0.02. Again, 
the stability of the solutions can be found by an examination of a linearised analysis of 
the growth rates of small perturbations to the equilibrium solutions. For K v < 0, the 
parameter a in equations (3.4a-c) is positive and from Segel and Stuart (l 96^), we deduce 
that the zero amplitude solution, I, loses stability at e = a corresponding to B2, so that a 
hexagonal pattern corresponding to either Ilia or IVa will be generated as e is increased 
further. This mode however is destabilised at A2. 

This loss in stability is due to the effects of the mixed mode V, which is unstable for 
all values of e, and originates as a bifurcation from the roll mode lib at Al. Beyond Al, 
lib becomes stable, and thus as the Rayleigh number is further increased, the hexagonal 
cells are replaced by the rolls corresponding to this solution. If the Rayleigh number is 
subsequently decreased, hysteresis can occur with the roll pattern persisting until Al at 
which point the hexagonal cells are re-established. 

For K p > 0, a similar scenario occurs. However, the parameter a is, in this case 
negative and thus instead of hexagonal cells corresponding to the equilibrium solutions 
Ilia or IVa being established at B2, the cell patterns Illb or IVb are initially generated. 
In addition the mixed mode bifurcates from rolls of type Ila, and it is this mode which 
becomes the stable solution beyond Al, with the roll pattern established at A2. 

In addition, it should be noted that the roll mode solutions cannot exist until e = —a, 
and for all values of K p , a is positive. This implies that the initial destabilisation of the 
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conduction mode, I, and transition to the hexagonal cells, occurs for values of e in excess 
of e = 0, the critical value of the uncontrolled problem. However, the hexagonal type 
solutions III, and IV , become viable solutions for t < 0 for all K p ^ 0, and thus if the 
Rayleigh number were decreased, hysteresis effects occur for subcritical values of e, that 
is, the closed cell pattern persists for negative values of this parameter. Figures 4 and 5 
describe how the location of the points B 1, B2 and Al, >12, vary with K p and illustrates 
the points discussed above. 

The stability of the various modes for K v > 0 and K p < 0 are outlined below. 


K p < 0 

K p > 0 

I 

Stable until 

e > —a 

Stable until 

e > —a 

Ila 

lib 

Unstable for all e 
Stable for 

e > Al 

Stable for 
Unstable for all t 

t > Al 

Ilia, IVa 
Illb, IVb 

Stable for 
Unstable for all e 

Bl < e < A2 

Unstable for all t 
Stable for 

Bl < e < A2 

V 

Unstable for all e 


Unstable for all e 



Figure 3, shows how the bifurcation picture develops for larger values of | K p |. As this 
parameter is increase, B 1 — ► — oo and B2 — ► -foo, and thus the no-motion state remains 
stable for longer. In addition, subsequent destabilisation of the hexagonal pattern is also 
delayed, that us, the range of e over which solutions III or IV are stable increases with 
K v 

We next consider the effects of the nonlinear controller on the bifurcation diagram, 
with the proportional controller, K p , set to zero. 

In figure 6, we show the bifurcation diagrams for small values of the nonlinear con- 
troller, K n = ±0.02. As in the previous case, stability of the various solutions depends on 
the sign of a in equations (3.4a). For, K n > 0, a is negative and vice versa. Thus the sta- 
ble/unstable modes are as previously discussed and outlined in Segel (1965). Again, after 
the loss of stability of the conduction mode I, an hexagonal cell pattern corresponding to 
either solution III or IV , is initially established. Due to the mixed mode V, these modes 
loses stability at A 2, and transition to the stable roll pattern II occurs as e is increased. 
If the Rayleigh number is then decreased hysteresis occurs, with the hexagonal mode not 
being re-established until Al. 

One important feature to be considered in this case is that both the roll and hexagonal 
modes become viable solutions for e < 0. Thus the initial transition from the no-motion 
state to the closed cells of solutions III and IV, occurs as a subcritical bifurcation. In 
figure 8 where we plot B\ and B2 , the thresholds above which solutions III, IV and II 
may exist respectively, against K n . As can be seen for all K n , these threshold values are 
negative, and thus transition from the no-motion state occurs for subcritical values of 
e. In addition, we see from figure 9, that Al and A2, the positions where the roll mode 
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becomes stable and destabilisation of the hexagonal cells occur respectively, occur for 
positive e. 

Finally let us look at the case where both controllers are in operation. Figure 10, 
illustrates the change in behaviour of the points 51 and 52, which represent the values 
of t above which solutions III, IV and II exist respectively. We recall that for K n 0, 
for all K p , 51 < 0 and 52 > 0, and for K p = 0, for any value of K n , 51, 52 were both 
negative. An important feature shown in figure 10 is that this behaviour can be modified 
such that there now exists a range of values of K p for a given K n , for which 51 > 0, 
that is to say the hexagonal cell pattern is not a viable solution until above e = 0, the 
critical value of the uncontrolled problem. In this case K n was varied between 0.0, and 
0.2 . Negative values of this parameter give similar plots, reflected in the line K p = 0.0. 

In addition, by careful choice of K p for a given K n , it is possible to create scenarios 
with, 

51 <0, 52 < 0, 51 <0, 52 > 0, 51 >0, 52 > 0. 

Since 51 < 52 for all K p , K n , it is not possible to generate the case where, 51 > 0, 52 < 
0. 

The bifurcation diagrams for these different cases are shown in figure 11. Here, solid 
lines represent stable solutions to the problem, whilst the dashed lines are unstable. In 
addition, we note that for the values in the parameters considered, except for certain 
cases, the solution corresponding to the mixed mode does not show up in these diagrams. 
This is due to the fact that the bifurcation from which this mode originates occurs for 
large values of e, out of the range illustrated here. However, behaviour will be similar 
to that discussed above, with the mixed mode originating as a bifurcation from the roll 
modes, at e = A\. Above this value, the roll mode solution will be stable, and the 
hexagonal cell structure will subsequently be destabilised due to interaction with the 
mixed modes, at A2. For larger e, either Ila or lib provides a stable equilibrium solution 
to the equations of motion. 

For K p = -0.55, K n = +0.05, we see that the hexagonal cell pattern becomes viable 
for e < 51 < 0, whilst the roll modes cannot occur until e > 52 > 0. This is due to the 
fact that the linear controller dominates over the nonlinear controller K n , and thus we 
see behaviour typical of the K n = 0, K p ^ 0 discussed previously. 

As K p is increased 51 moves through the origin, and solutions III and IV occur 
for positive value of e, whilst 52, moves towards t — 0. This behaviour is due to the 
relative dependence of the parameters ct and c on K p and K n . For this parameter range, 
some cancellation of the effects of each controller on a and a is occurring leading to 
supercritical behaviour, that is to say that the destabilisation of the no-motion state I, 
occurs for values of e in excess of the critical value of the uncontrolled problem. These 
points eventually coalesce and this case, being of some interest, is discussed later. 

For further increases in K p , 51 again becomes negative and, for even larger K p the 
roll solutions II also become viable for e < 0, and thus the initial bifurcation form 
the conductive state to the hexagonal cell pattern is subcritical in nature. Within this 
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parameter range, the nonlinear controller dominates the regulation of the lower surface 
temperature. B2 again becomes positive for larger values of K p , with B 1 remaining 
negative in value. Here the relative effects of the controllers are such that the linear 
controller again dominates. 

The stability characteristics for each of these cases are as previously discussed. The 
no-motion state loses stability at e = -a, corresponding to point B2, and the hexagonal 
cell pattern, forms as the stable solution. Due to the influence of the mixed mode V, 
which results as a bifurcation from the roll modes, this pattern is then destabilised at 
A2, and rolls are formed. Subsequent decreases in the Rayleigh number lead to hysteresis 
effects, with the roll modes persisting until Al, at which point the hexagonal cells are 
re-established. 

As mentioned previously, points B\ and B2 can coalesce. From the form of the 
equilibrium solutions given in (3.4a-e), we can see that this occurs when a = 0. Due to 

the dependence of a on K p and K n , given by, (2.26), a relationship between the control 
parameters maybe written down, 


K, = 

d\ 


(3.7) 


where a, and a 2 are constants. In addition at this point, the equilibrium solutions reduce 
to, 


I Conduction : x 0 = y 0 = z 0 = 0, ( 3 _g a ) 

IIa,b Rolls : x 0 = y 0 = 0, z 0 = ± 2 , (3.8b) 

Ilia, b Hexagons : x 0 = y 0 = z 0 , z 0 = T (2T)~ l ^4 (e + a) T, (3.8c) 

IVa,b Hexagons : z 0 = Vo = ~z 0 , z 0 = zf (2T)~ l ^/4(e + a)7\ (3.8d) 

V Mixed : z 0 = 0, z 0 = y 0 = ±y/R(e + a). (3.8e) 


A typical bifurcation diagram for K n and K p satisfying (3.7), is shown in figure 12. 
Here, solutions corresponding to IIa,b are shown as dotted lines, whilst the hexagonal 
cell patterns IIIa,b and IVa,b are solid. As can be seen, the picture is very similar to that 
for the equilibrium solutions for * 0 with no control, (figure 1), except that the bifurcation 
from the no-motion state occurs for e = -a. The roll and hexagonal modes become 
possible solutions to the equations of motion at this point, as indeed does the mixed 
mode, which corresponds to a solution with no-motion in the ^-direction. Beyond e = a. 
the conduction mode, loses stability, and the hexagonal cell pattern is formed. Since the 
mixed mode V originates as a bifurcation from the roll solutions II at the same point, 
e = a, and no interaction between this mode and modes III or IV is possible, there can 
be no subsequent transition from the hexagonal cells to rolls . In addition, if the Rayleigh 
number were decreased from above e = -a, no hysteresis effects are seen. 
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4 Further comments and conclusions 


We have given a description of the process by which a control strategy maybe used to 
delay or indeed speed up the onset of convection in a Boussinesq fluid. The amplitude 
equations derived describe the effects of both linear and nonlinear proportional controllers 
on the flow regime. 

Since the controller has no preferred direction, a subcritical bifurcation to the hexag- 
onal cell pattern remains viable. We considered three cases in detail. Firstly, with only 
the linear controller in operation, we found that the onset of convective instabilities was 
delayed, and that the hexagonal cell pattern persists for longer as the control parameter 
K p is increased. If the linear controller is switched off, and thus the lower surface tem- 
perature is only regulated by the nonlinear controller, the transition from the no-motion 
state to a hexagonal cell pattern occurs for subcritical values of the Rayleigh number. 
The most interesting case was with both controllers in operation. Here, a range of dif- 
ferent behaviours is possible. In particular the exists a parameter range for which the 
hexagonal cell patterns do not become viable solutions to the model for values of e greater 
than zero. Thus, not only is the onset of convection delayed significantly, but hysteresis 
effects would lead to the re-establishment of the no-motion state for values of e in excess 
of the critical value for the uncontrolled problem. In addition, for a given value of the 
parameter K n , we have shown that there exists a K p such that the onset of convection is 
supercritical, and the solutions corresponding to the roll, hexagonal and mixed modes all 
bifurcate from the same point, thus the hexagonal cell pattern remains the stable solution 
for Rayleigh numbers in excess of this point. 

To conclude, careful choice of the various parameters involved, can either delay the 
onset of convection, or produce transition for subcritical values of the Rayleigh number. 
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